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O 1 ABSTRACT 

Context. Earth-like planets have viscoelastic mantles, whereas giant planets may have viscoelastic cores. The tidal dissipation of such 
\Q • solid regions, gravitationally perturbed by a companion body, highly depends on their rheology and on the tidal frequency. Therefore, 

modelling tidal interactions presents a high interest to provide constraints on planets' properties and to understand their history and 
i— i , their evolution, in our Solar Sytem or in exoplanetary systems. 

Qi^ , Aims. The purpose of this paper is to examine the equilibrium tide in the anelastic parts of a planet whatever the rheology, taking into 

■ account the presence of a fluid envelope of constant density. We show how to obtain the different Love numbers that describe its tidal 
deformation. Thus, we discuss how the tidal dissipation in solid parts depends on the planet's internal structure and rheology. Finally, 

i-C ' we show how the results may be implemented to describe the dynamical evolution of planetary systems. 

Methods. We expand in Fourier series the tidal potential, exerted by a point mass companion, and express the dynamical equations 
JL 1 in the orbital reference frame. The results are cast in the form of a complex disturbing function, which may be implemented directly 

r* l ' in the equations governing the dynamical evolution of the system. 

, Results. The first manifestation of the tide is to distort the shape of the planet adiabatically along the line of centers. Then, the 

response potential of the body to the tidal potential defines the complex Love numbers whose real part corresponds to the purely 
i i ■ adiabatic elastic deformation, while its imaginary part accounts for dissipation. The tidal kinetic energy is dissipated into heat through 

anelastic friction, which is modeled here by the imaginary part of the complex shear modulus. This dissipation is responsible for 
7-H | the imaginary part of the disturbing function, which is implemented in the dynamical evolution equations, from which we derive the 

^ , characteristic evolution times. 

Conclusions. The rate at which the system evolves depends on the physical properties of tidal dissipation, and specifically on how 
\Q ' the shear modulus varies with tidal frequency, on the radius and also the rheological properties of the solid core. The quantification 

of the tidal dissipation in solid cores of giant planets reveals a possible high dissipation which may compete with dissipation in fluid 
t-H . layers. 

Key words. Planet-Star interactions - Planets and satellites: general - Planets and satellites: physical evolution - Planets and 

■ satellites: dynamical evolution and stability 
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1 . Introduction and general context the orbital motion) or the companion tends to spiral into the par- 
ent body. Indeed, by converting kinetic energy into heat through 

Since 1995 a large number of extrasolar planets have been dis- internal friction, tidal interactions modify the orbital and rota- 

H ■ covered, which display a wide diversity of physical parameters tional properties of the components of the considered system, 

- - (Santos & et al., 2007). Quite naturally the question arose of their and thus me ir structure through internal heating. This mecha- 

habitability, i.e. whether they could allow the development of nism depends sensitively on the internal structure and dynamics 

life. Determining factors are the presence of liquid water and of f tne perturbed body 

a protective magnetic field, properties which are closely linked Recent smdies haye been carried QUt Qn Mal effectfJ in fluid bod . 

to the values of the rotational and orbital elements of the plane- ies nke stars and envel of iant lanets (Q iMe & Lin 20 04- 

tary systems. And these elements strongly depend on the action 2007; Qgilvie 2009; Remus, Mathis & Zahn 2012). However 

of tides. Once a planetary system is formed in a turbulent accre- ±e planetary solid regionSi such as mantles of Earth . lik e plan- 

tion disk, its fate is determined by the initial conditions and the etfj Qr rQck CQres of iant lanetS) if preS ent (e.g. Guillot 1999, 

mass ratio between planet and hosting star. Through tidal inter- Gaulme gt aJ 2Q1 ^ contri bute likewise to tidal dissipa- 

actton between components, the system evolves either to a stable tion The firfjt mdy of a tidally deformed elastic body was done 

state of minimum energy (where all spins are aligned, the orbits b Lord KeMn ( { g63) who Hed k tQ an incompressible ho _ 

are circular and the rotation of each body is synchronized with mogeneous Earth . Further developments were made by Love 

(1911), who introduced a set of dimensionless numbers, the so- 

Send offprint requests to: F. Remus called Love numbers, to quantify the tidal perturbation. More 
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recently Greff-Lefftz (2005) generalized these results in the case 
of a spheroidal rotating Earth. In the meanwhile Dermott (1979) 
considered a two-layer model and he studied the impact of a 
tidally deformed static fluid shell on the deformation of an elas- 
tic solid core. 

If the body is not perfectly elastic, i.e. if its internal structure 
is anelastic, the tidal deformation presents a lag with respect to 
the tension exerted by the tidal force, and causes the dissipation 
of kinetic energy. Several recent studies addressed this prob- 
lem of tidal dissipation using linear viscoelastic models. Peale 
& Cassen (1978) evaluated the tidal dissipation in the Moon 
considering various models of internal structure. Tobie et al. 
(2005) applied the Maxwell rheological model to evaluate the 
dissipation in Titan and Europa; Ross & Schubert (1986) com- 
pared three different linear models of viscoelasticity (Kelvin- 
Voigt, Maxwell, Standard Anelastic Solid), to which Pfenning 
et al. (2009) added the Burgers body. All these studies suscitate 
our interest in the tidal dissipation resulting from the anelastic 
deformation of the solid parts of a planet when perturbed by a 
companion. 

We shall study here the tidal dissipation in a planet which pos- 
sesses an anelastic core, made of a mix of ice and rock, sur- 
rounded by a fluid envelope, such as an ocean. The planet is 
part of a binary system where what we call the companion (or 
perturber) may be either the hosting star or a satellite of the 
planet. Due to the tide exerted by the companion, the core of 
the two-layer planet is deformed elastically, but because of the 
anelasticity of the material composing the core, this deformation 
is accompanied by a viscous dissipation that we evaluate what- 
ever the rheology. As an illustration, the results will be given 
for a Maxwell body. We then compare the value of the tidal dis- 
sipation in presence of a fluid envelope with that achieved by 
the fully solid planet, and we examine the dependence of the re- 
sults on the relative sizes of the core and the planet, the relative 
densities, and the viscoelastic parameters. In the last section, we 
establish the equations governing the dynamical evolution of the 
system, from which we deduce the caracteristic times of circu- 
larization, synchronization and spin alignments. 



2. Elastic deformations of a solid body under tidal 
perturbation 

2. 1 . The system 

We consider a planet A of mass Ma, consisting of a rocky (or icy) 
core and a fluid envelope, rotating at the angular velocity Q and 
tidally perturbed by a second body of mass Mb, assumed to be 
ponctual, moving around A on a Keplerian orbit, of semi-major 
axis a and eccentricity e, at the mean motion n. We locate any 
point M in space by its usual spherical coordinates (r, 9, ip) in 
a spin equatorial reference frame %e '■ {A,Xe, Ye,Ze) centered 
on body A and whose axis (A, Ze) has the direction of the rota- 
tion axis of A. The corresponding configuration is illustrated on 
Figure 1. 

In this section, we first assume that planet A has no fluid layer 
and its internal structure is supposed to be perfectly elastic. We 
will then denote by p its density and R its mean radius. 



2.2. The tidal potential 

The planet is submitted to a tidal force, exerted by the perturber 
B, which derives from a perturbing time dependent potential 




Fig. 1. Spherical coordinates system attached to the equatorial 
reference frame 'Re ■ {A,Xe,Ye,Ze} associated to body A. A 
point M is located by r s (r,Q,ip); the point mass body B by 
*b = (r B , B , <Pb)- 



U(r, t). Following Zahn (1966a-b) and generalizing it by using 
Kaula (1962), Lambeck (1980), Yoder (1995-1997) and Mathis 
& Le Poncin (2009) (hereafter MLP09) in the present case of a 
close binary system where spins are not aligned, the components 
are not synchronized with the orbital motion and where the or- 
bit is not circular, we expand the tidal potential U in spherical 
harmonics Y'"(6, if) in 'Re- 

Before we proceed, we need to define the Euler angles 
that link the spin equatorial frame He ■ {A,Xe,Ye,Ze} of 
the central body A, on one hand, and the orbital frame 
%o '■ {A, Xo, Yo, Zo), on the other hand, to the quasi-inertial 
frame %r : {A, Xr, Yr, Zr) whose axis Zr has the direction of 
the total angular momentum of the whole system. 

We need the three following Euler angles to locate the orbital 
reference frame %o with respect to 

- I, the inclination of the orbital plane of B; 

- of, the argument of the orbit pericenter; 

- £2*, the longitude of the orbit ascending node. 

The equatorial reference frame %e is defined by three other Euler 
angles with respect to Kr: 

- s, the obliquity of the rotation axis of A; 

- 0*, the mean sideral angle defined by Q. = d0*/dr; 

- <p*, the general precession angle. 

Refer to Figure 2 for an illustration of the relative position of 
these three reference frames and the associated angles. For con- 
venience, all the following developpements will be done in the 
spin equatorial frame %e of A (as it has been done in MLP09). 

All following results are derived from the Kaula's transform 
(Kaula 1962), used to explicitly express the whole generic multi- 
pole expansion in spherical harmonics of the perturbing potential 
U in terms of the Keplerian orbital elements of B in the equato- 
rial A-frame: 
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B(a„,e B ) 




periastron 



Fig. 2. Inertial reference CRr), orbital (7?o), and equatorial (He) 
rotating frames, and associated Eulers angles of orientation. 



B j=-l p=0 qeZ 



2i + i(i-\j\y. 

4n + 



xd\ m {e)F Up (I)G, 4 ,, q (e)e^ 



where 9b and tpB are respectively the colatitude and the longitude 
of the point mass perturber B, and where the phase argument is 
given by: 

^l^n,j,p,q(t) = CTl, nhPt q(n, Q.) t + Tf m jpq (<j/,Q*, f) . (2) 

We have defined here the tidal frequency: 

o~i,m, P ,q{n, Q) = (I - 2p + q)n- m£] , (3) 
and the phase T; mjjP 9 : 

Ti, m ,j, P , q = (/ - 2p)aj* + j(Q.* - f) + (/ - m) |. (4) 

We study here binary systems close enough for the tidal interac- 
tion to play a role, but we also consider that the companion is far 
(or small) enough to be treated as a point mass. We then are al- 
lowed to assume the quadrupolar approximation, where we only 
keep the first mode of the potential, I = 2: 



!2 2 2 
z zzz 
m=-2 7=-2 p=0 qeZ 



where 



U m , jtM {r)P^{co^G)e'^ 



^>2,m,j,p,q(lfi, t) = mtf + y ¥2,mj,p,q(t) 



,<4>2,„,j, M (¥>.') I 



(5) 



(6) 



The functions U m j iPiq (r, 9) may be expressed in terms of the 
Keplerian elements (the semi-major axis a of the orbit, its eccen- 
tricity e and its inclination T) and the obliquity e of the rotation 
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Table 1. Values of the obliquity function d 2 m (e) in the case 
where j > m obtained from Eq. (8), (from MLP09). 



axis of A, as: 

U m ,j, P , q (j) = (-1)" 



(2-m)!(2-|j|)! 
^(2+m)!(2 + |;|)! 



d\ m (e)F %k ,(l)G 2 , m (e)\ r 2 , (7) 

where Q is the gravitational constant. 

The obliquity function d 2 m {y) is defined, for j > m, by: 



4m(r) = (-ir m 
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where P ( (x) are the Jacobi polynomials (cf. MLP09). The val- 
ues of these functions, for indices j < m, are deduced from 
- (— l) 2- ^d^ m (y) or from their symmetry proper- 
: {—\y~'"d 2 . (y) = d 2 (-y); moreover, we have: 

2 



ties: d 2 m (y) 

= 5j jn . Values are given in Table 1. 
We also define, the inclination function F2j lP (I): 
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with the symmetry property: 



F2,-j,p(D 



(-1) 



F 2,j,p (J) ■ 



(10) 



(2 + ;)!_ 

Values are given in Table 2. 

The eccentricity functions G2, P , q {e) are polynomial functions 
having e q for argument (see Kaula, 1962). Their values for the 
usual sets {2, p, q] are given in Table 3, knowing that in the case 
of weakly eccentric orbits, the summation over a small number 
of values of q is sufficient (q e [ - 2, 2J). In the following, let us 
denote by I = [ - 2, 2] x [ - 2, 2] x [0, 2] x Z the set in which 
the quadruple \m, j, p, q] takes its values. 

If we simplify the expansion of the potential in the case where 
spins are aligned and perpendicular to the orbital plan, where 
obliquity s and orbital inclination I are zero, Eq. (5) reduces to 
the expression of the potential given by Zahn (1977). 
The tidal force induces a displacement of each particule consti- 
tuting the planet, thus causing some deformations. In particular, 
the core's surface is deformed as described by the Love theory 
(Love 1911). 
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Table 2. Values of the inclination function F2j, p (I). Values for 
indices j < can be deduced from Eq. ( 10), (from MLP09). 
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Table 3. Values of the eccentricity function G2, p , q (e), (from 
MLP09). 



2.3. Dynamical equations for a solid body and their boundary 
conditions 

To describe the internal evolution of the main component A sub- 
mitted to the perturbations induced by the tidal potential pre- 
sented above, we use the Eulerian formalism (Dahlen et al. 
1999). The system of equations, needed to follow the motion of 
a particule, is composed by the Eulerian momentum (11a) and 
mass (lib) conservation laws, and the Poisson equation (11c) 
satisfied by the potential <£ of self-gravitation: 



dt 2 



V-o-+pV(<D + f/) 



-i- + V • In— = . 

dt \ H dt J 

V 2 <D = -AitQp , 



(11a) 

(lib) 
(He) 



where s designates the displacement vector and & is the stress 
tensor. We complete this system with the constitutive equation to 
link the stress exerted on the body to the resulting deformation. 
Assuming that tidal effect corresponds to a traction applied on 
the body, with no rotational contribution, the deformation tensor 
reduces to the strain tensor e: 



where h designates the transposed tensor of h. We then get a 
relation linking the stress tensor & to the strain tensor e that ac- 
counts for the rheology of the body, and that we represent by a 
function f r h- 

B- = T A {e). (He) 

To solve this system (11), we need to apply boundary condi- 
tions to the five previous equations, assuming that there is no 
displacement (12a) neither attraction (12b) at the center of mass 
r — 0, the gravitational potential has to be continuous (12c) and 
the Lagrangian traction has to vanish (12d) at the surface r - R: 



s| r= o=0, 

(O + W\r=0 

[<»]£ = 



50 A ^ 

— +\nQps r 
or 



= 0. 



(12a) 
(12b) 

(12c) 
(12d) 



2.4. Linearization of the system 

Assuming that tidal effects, and thus the resulting elastic de- 
formation, are small amplitude perturbations to the hydrostatic 
equilbrium, we are allowed to linearize the system (11) and its 
boundary conditions (12). To do so, we expand a scalar quantity 
Xas: 

X (r, 6, ip, t) = X (r) + X' (r, 0, <p,t) ; (13) 

Xq designates the spherically symmetrical profile of X, and X' 
represents the perturbation due to the tidal potential. The dis- 
placement s and the tidal potential U are also considered as per- 
turbations. Thus, correct to first order in ||s||, we obtain the fol- 
lowing form of system (11): 

dh 
Po dfi 



where we made use of the Hooke's law (14d), which is a lin- 
ear constitutive law that governs elastic materials as long as the 
load does not exceed the material's elastic limit, in the case of 
an isotropic material (i.e.: whose properties are independent of 
direction in space). It means that strain is directly proportional 
to stress, through the bulk modulus K and the shear modulus /j. 
The reference state is drawn from an up-to-date planetary struc- 
ture model. It is governed by the following Poisson equation and 
the static momentum equation: 



V ■ 0- + po V (O' + U) + p V<E>o , 


(14a) 


0. 


(14b) 


-4nQp' , 


(14c) 






(tf-|ultr(i)l + 2/ii, 


(14d) 



V 2 <D = -AnQpo , 
VP =PoVcDo, 



(15a) 
(15b) 



= \ [Vs + (Vs) T ] , 



(Hd) 



where we made use of the following convention for the gravity: 
go = V(D . 

2.5. Analytical solutions for an homogeneous incompressible 
body 

To solve the linear system (14), we expand all scalar quan- 
tities in spherical harmonics Yf(6,(p). Moreover, as all vecto- 
rial quantities that intervene in Eqs. (14a-14b) are poloidal, we 
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may expand them in the basis of vectorial spherical harmonics 
[R™(0, <p), Sf(0, 9?)], where R refers to the radial part and S to the 
spheroidal part of a given vector (Rieutord 1987, Mathis & Zahn 
2005): 



V(/,m)£NxI-/ ) /], 



R™(G,<p) = Y?(G,<p)e r , 

s';\0,<p) = v 5 [y?<{0,<p)] , 



(16a) 
(16b) 



where Vs designates the horizontal gradient: 



d- 1 a- 
89 sin 6 dip 



(17) 



We introduce six radial functions y"[ 6 j(r) (Takeuchi & Saito 
1972) to expand all quantities in spherical harmonics, at the 
quadrupolar approximation (Z = 2): 



- the displacement: 



s(r,e,<p,t)= Yi [yTtoH?^) 



(m,j,p,q)€l 



+y%(r) S%(6, if)] e^W-W , (18a) 



- the total potential: 



lutions, considering the expansion (18): Vm e [ - 2, 2], 



k 2 /8 



r 2 R gs \3 



2/7 



+ -nQp L -^- \-R--r 



ki /8 n2 1 



■p(l+k 2 ) 



- V 3 



- V 4 



- v 5 



n 



5 , 



Um,j,p,q(r) ■ 



(r) = X ^ -j^- (# 2 - r 2 ) C/ mJ>Pj? (r) , 

W = J](l+fe)C/m,;; M W, 



2(1 +/c 2 ) 



j,p,q 



r 

— I 



^ T m,j,p,q(y) ■ 



(19a) 



(19b) 

(19c) 

(19d) 
(19e) 



(19f) 



where we have introduced the acceleration of gravity at the sur- 
face g s and the second-order Love number k 2 . The latter com- 
pares the perturbed part <&'(R) of the self-gravitational potential 
at the surface of a fully-solid planet of mean radius R, deformed 
by tidal force, with the tidal perturbing potential U(R): 



(U+®')(r,e,<p,t) = Y y';(r)Y'l\e,ip)e fV2 ^< (f) , (18b) 

(m,j,p,q)£l 



- the Lagrangian traction: 



def &(R) 

" U(R) 



(20) 



The expression of k 2 is established in Sect. 3, for an ocean-free 
planet (Eq. 61) or a two-layer planet (Eq. 63). 



e r -&(r,0,tp,t)= J] [y2(r)K2(e,<p) 

(m,j,p,q)£l 

+y'Z(r)S'2(8,<p)]e i ^ W> , (i 8c ) 



- the Lagrangian attraction (introduced to express the continu- 
ity of the gradient of the potential): 



Vm e I - 2, 2], y»(r) = [)£(r)] - 4^p ^(r) . (18d) 



The linear system governing the radial functions yj[ 6 j(r) is 
given in Appendix. In the case of an incompressible (K — > +oo) 
and homogeneous body (p, po = csfX the system (14) con- 
strained by boundary conditions (12) admits the following so- 



3. Modified elastic tidal theory in presence of a fluid 
envelope 

We now assume that planet A is not entirely solid, but has a fluid 
envelope. We follow the method proposed by Dermott (1979) to 
evaluate how the anelastic dissipation is modified by the pres- 
ence of a fluid layer surrounding the solid region. The first step 
consists in determining the behaviour of the elastic response in 
this configuration. We denote by R c (resp. R p ) the mean radius 
of the solid core (resp. of the whole planet, including the height 
of the fluid layer); p c and p designate the density of the core 
and the ocean, both assumed to be uniform, as a first step. More 
generally, all quantities will be written with a "c" subscript when 
evaluated at the core boundary, and with a "p" subscript if taken 
at the surface of the planet. The evolution of the system is de- 
scribed in the orbital frame 'Ro ■ {A,Xo,Yo,Zo} centered on A 
and comoving with the perturber B. We will use polar coordi- 
nates (r, 0) to locate a point P, where r is the distance to the 
center of A, and © is the angle formed by the radial vector and 
the line of centers. 
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Fig. 3. Left : tidal displacement s. Middle: equatorial slice of s. Right: meridional slice of s. The orange arrow indicates the direction 
of the perturber B, the red one corresponds to the rotation axis of A. The two slices are planes of symmetry. 



3.1. Vertical deformation at the boundary of the core 
In %o, the tidal potential takes the form (Dermott 1979) 

2 

U(r) = -ar)g(r)P 2 (cos&) = -( cgc — P 2 (cos0) , (21) 



where we have introduced the tidal height 



M, 



and the gravity 



g(r) = 



— (-f 

M(r) W 

GM{r) 



(22) 



(23) 



M(r) being the fraction of mass of the planet inside the radius r. 
The expression of the tidal potential in the rotating frame of 
B (Eq. 21) is linked to its expression in the equatorial inertial 
frame (Eq. 5), through the Kaula's transform (Eq. 1). Indeed, the 
Legendre polynomia summation formula 



P 2 (cos 0) = fk 



1 

m=0 



(2-m)! 
(2 + m)! 



(2 -£„,„,) P'^cos^e 



x (cos Gb) e 



tlllif 



(24) 



involves the term P™ (cos B ) e"' m ^ in Eq. (21) that has to be 
transformed following Eq. (1-2) to obtain Eq. (5). 

In this section we are interested in the mofication of Love num- 
bers due to the presence of a fluid envelope on top of the solid 
core. Thus, we will focus on the deformations of the core's 
surface, and particularly on the vertical displacements. Love 
(1911) proved that tidal deformations could be described by the 
same harmonic function than the tidal potential which causes it. 
Therefore the equations of the core and planet boundaries are 
respectively of the form 

r c = R c + s r (R c ) = R c [1 + S 2 P 2 (cos ©)] , (25a) 
r p =R p [l + T 2 P 2 (cos®)] . (25b) 

Thus, s r {R c ) = R c S 2 P 2 (cos 0) represents the radial displace- 
ment at the core's boundary corresponding to the vertical tidal 
deformation of amplitude S 2 P 2 (cos®). In 1909, Love defined 



the number h 2 , as the ratio between the amplitude of the verti- 
cal displacement at the surface of the planet and the equilibrium 
tidal height (disturbing potential divided by undisturbed surface 
gravity, both taken at the surface of the core) in the case of a 
fully-solid planet. Solving the whole system of equations, he de- 
termined its expression as 



hi = 



def S r (R c ) R C S 2 5 1 



U(R c )/g c 



Zc 2 1+/2* 



(26) 



where fi is called the effective rigidity, in the sense that it evalu- 
ates the relative importance of elastic and gravitational forces: 



f = 



19// 



(27) 



In presence of the fluid envelope, the ratio between the amplitude 
of the tidal surface vertical displacement and the tidal height will 
be modulated by a multiplicative factor F, due to the additional 
loading exerted by the tidally-deformed fluid layer. We may then 
introduce a new notation h^ for the modified Love number in 
presence of a fluid envelope: 



Sr(Rc) 

U(R c )/g c 



Fh 2 = Fx 



5/2 



1 



(28) 



We have now to express this factor in function of the parameters 
of the system. To do so, we have to list all the forces acting on 
the surface of the core. Before carrying out the study of these 
forces, let us introduce a specific notation. All physical quanti- 
ties X(r) will be separated in two terms: the first corresponds to 
the constant part that does not depend on where the quantity is 
calculated; the second one (called the "effective deforming" con- 
tribution and denoted X'(r)) is a term proportional to the spheri- 
cal surface harmonic P 2 (see 24). 

3.2. Gravitational forces acting on the surface of the core 

The planet is not only subjected to the direct action of the tidal 
potential U, but also to the self-gravitational potential O per- 
turbed by the first. In calculating the latter, we have to consider 
both contributions of the solid core and the fluid envelope, <t> c 
and <t> respectively. 

At any point r of the core, O^r) corresponds to the internal po- 
tential created by the core : 



3Rt - r 2 



+ -SS 2 P 2 



(29) 
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At the same point r, 0„(r) is the internal potential created by the The perturbations of pressure P'(r, 0) obey the relation of the 



fluid shell of density p D and of lower and upper surface bound- 
aries r c and r p respectively: 



O c (r) 



Po gc_ 
Pc Rc 



(3R 2 p -r 2 3 



3R 2 -r 2 3 , 

— ^ + -=r L S 2 P 1 

2 5 



(30) 



Therefore V(r) = U(r) + <£ c (r) + <D (r) has the following ex- 
pression, at any point r inside the core: 



V(r) = -- & 

2 R r 



-r + 3Rt 



K>2\ 



Pc Pc R 2 c, 

where the effective deforming potential is expressed by 
V(r) 



V'(r), Ola) 



-Zr 2 Pi 



Z being a constant that depends on the characteristics of the 
planet: 

Z=ii k. + lPl(T 2 -S 2 )+ -S 2 . (31c) 
R c [R c 5 Pc 5 J 

We then obtain its expression, correct to first order in S 2 P 2 or 
T 2 P 2 , at any point r c = r c e r of the surface of the core: 



V(r c ) 
where 



ScRc 



1 _ ISA + -El 
2p c ) 2p c R 2 C 



V'(r c ) = -ZcR 2 c P 2 , 



V'(r c ) 



(32a) 



(32b) 



Z c being a constant that depends on the characteristics of the 
planet: 



Z c = — 

Rr 



k + \ P -{T 2 -S 2 )--S 2 
R c 5 p c 5 



(32c) 



Chree (1896) showed that the deformation of the core's surface 
under the gravitational forces (that derive from the effective de- 
forming potential V) is the same as that which would result 
from the outward normal traction f x N applied at its mean sur- 
face r = R C : 

f^(R c )=PcZR 2 c P 2 . (33) 

3.3. Total effective normal traction acting on the surface of 
the core 

The mean surface of the core is subjected to two additional 
forces induced by both the loading of the core and the loading of 
the ocean, tidally deformed. 

First, the pressure due to the differential overloading of the de- 
formed elastoviscous matter on the mean surface of radius R c is 
given by the product of the local gravity g c , the density of the 
core p c and the solid tidal height R C S 2 : 



ft N (Rc) = -PcgcRcS 2 P 2 



(34) 



We also have to take into account the oceanic hydrostatic pres- 
sure. Following Zahn (1966) and Remus et al. (2012), we ex- 
press all scalar quantities X(r, 0) as the sum of their spherically 
symmetrical profile X°(r) and their perturbation X'(r, 0) due to 
the tidal potential U(r, 0) oc P2(cos 0): 

X(r, 0) = X°(r) + X'(r, 0) = X°(r) + X{r) P 2 (cos 0) . (35) 



hydrostatic equilibrium which is of the following form, correct 
to first order in P2(cos 0): 



VP' = Po W + p' D vv° . 
Therefore, the 0-projection of (36) leads to 
P'(r,®)=p V'(r,®). 



(36) 



(37) 



Finally, since only the variable part of the pressure, i.e. P', con- 
tributes to the normal effective traction fI N that acts on the mean 
surface of core, this latter takes the following expression: 



€ N (R C ) = P'(r c ) = Po V'(r c ) = - Po Z c R 2 c P 2 



(38) 



The sum of these three forces, represented on Fig. 4, 
corresponds to the total normal effective traction 

(31b) f T HRc)=f( N (Rc) + f 2 N (Rc)+fi N (Rc) that deforms the 
mean surface of the core. Using Eqs. (33-34-38), we get: 



f N (R c ) = (pcZ-p Z c )R 2 cP 2 - Pc gc Rc S 2 P 2 , 



(39) 



where the expressions of Z and Z c are given by equations (31c) 
and (32c) respectively, so that 



f»(R c )=XP 2 (cos@), 



where we have denoted by X the following quantity 
X 



2 p , P» 

TPcgcRclt 

5 \ Pc 



Ik 

2Rc 



S 2+ 3 - P -(T 2 -S 2 ) 
2pc 



(40) 



(41) 




Fig. 4. Balance of forces that act on the mean surface of the core 
r — R c : f^ N (R c ) are the gravitational forces, f^ N (R c ) the loading 
of the solid tide and fj N (R c ) the hydrostatic pressure. 



3.4. Amplitude of the vertical deformation 

According to Melchior (1966), a deforming potential H 2 of sec- 
ond order produces a deformation at each point r c of the surface 
of the core which radial component takes the form 



{8R 2 .-3r 2 ) Pc ^ 2 



I9p 



(42) 
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Correct to first order in P 2 , recalling that r, 
reduces to: 
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R c (1 + S 2 P 2 ), this where 



(43) 



where p c 1l 2 is the deforming traction f T ' N (R c ) = XP 2 applied on 
the core. Furthermore, since the amplitude of the displacement 
is also given by e rr = S 2 P 2 (Eq. 25a), we have the following 
equality 

S 2 P 2 =-^-XP 2 . (44) 
19 p 

Therefore the relation (27) beetween p. and fi and the expression 
(41) of X enable us to relate the deformation of the surfaces of 
the core (S 2) and of the ocean (T 2 ) as 



P 



i-es. 



2R C 2p c 



(45) 



By definition, given in Eq. (28), the impedance F is of the form 

P=^(l +/ u)^S 2 . (46) 

J be- 



Since the surface of the planet is an equipotential, the total po- 
tential V takes a constant value at any point r p of the surface of 
the ocean: 

V(r p ) = U(r p ) + 4>(r p ) = est . (47) 
As V — V = est by definition, we get the simplier condition: 



V'(r p ) = t/(r p ) + 0'(r p ) = est . 



(48) 



At a point r of the ocean, <l> c (r) corresponds to the external po- 
tential created by the core: 



® c (r) = -g c Ri\- + - 



1 3R 2 C S 2 



(49) 



At the same point r, <t>„(r) is the internal potential created by the 
fluid shell of density p and of lower and upper surface bound- 
aries r c and r p respectively: 



Rc Pc 



(3 R 2 p -r 2 | 3 r 2 T 2 ^ 
2R\ 5 Rl 



2 p 1 1 3 R 2 S 2 



gcK- - + 7 

Pc 



Pi 



(50) 



Therefore V(r) = U(r) + <D c (r) + <t> (r) has the following ex- 
pression, at any point r inside the ocean: 



V(r) = -g c R z c 



1 -8l\ I _ Pf.— + -PjL^L 

Pc) r p c 2Rl IpcRl 



+ V'(r) , 

(51a) 



where the effective deforming potential is expressed by 

V'{r) = -g c R 2 c W{r)P 2 , (51b) 
W being a function of the distance r to the center of the planet: 



uif ^ & r 2 3 p„ r 2 3 
^ (r) ^^ + 5^^ r2 + 5 



(51c) 



We then obtain its expression at a point r p = r p e r of the surface 
of the planet: 



V(r p ) = 



1 _ Pi) ^£ + PlL^L 



Rn 



Pc Rc 



+ V'(r p ) , 



(52a) 



V'(r p ) = -g c W p P 2 , 



(52b) 



W p being a constant that depends on the planet characteristics: 



W„ 



5 

Rl 



Po 
Pc 



Rl 



2 


Po 


4 




SI 


~5 


Pc 


Rc 





T 2 . (52c) 



The condition (48) takes then the form: 



i = 4^i 11 



5 \R 



— I S 2 + - —aT 2 , 
5 p c 



where 



2 Po \R 



(53) 



(54) 



We may eliminate the variable T 2 thanks to Eq. (45): 



2 Po 7 
a 1 7 — 



5Pc ^ (a + l)(l-%) 



Po 3 p 



+ -— 1 +- — 



Pc 2pc 



Po 



3 IR C 



2 \R 



l- P -l\ 

Pc) 



(55) 



Injecting this relation into the expression (53) of £ C /Rc, and the 
resulting relation in the expression (46) of F, we finally get: 



(1 -£)(!+/» (l + £) 



(56) 



In the case of a shallow oceanic envelope (R p 7? c ), the height 
of the oceanic tide is then given by /? c (T2 - S 2 ) at the surface of 
the core. Using Eqs. (53-45), we obtain the classical expression 
of the height of oceanic tide 



Rc(T 2 -S 2 ) 



Pic 



(57) 



The height of the solid tidal displacement is given by R C S 2 . 
Using Eqs. (57-53-45), we obtain its classical expression: 



R C S 7 — 



which reduces to 



2 (p 
1 + p 



(58) 



(59) 



for an ocean-free planet (p a = 0), which corresponds to that 
given by Lord Kelvin (1863). Thus, recalling Eq. (46), we de- 
duce that for an oceanless planet, F is unity. 
Fig. 5 displays the value of F for three types of planets (i.e. 
Earth-, Jupiter- and Saturn-like planets), with a given core (of 
fixed size, mass and shear modulus) and a fluid shell of fixed 
density but variable depth, so that the size and mass of the whole 
planet varies also. The variation of F is represented in function 
of R c /Rp- the smaller this ratio, the higher the ocean depth. 
For the Earth, all parameters are well known (see Table 4). 
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The values of the ocean density for the Jupiter- and Saturn-like 
planets correspond to the ones we may deduce from the well 
known values of their global size and mass (see Table 5), and the 
much less constrained values of the core size and mass (see mod- 
els A of Table 6). Concerning the shear modulus, we used the 
value taken as reference by Henning et al. (2009) when studing 
the tidal heating of terrestrial exoplanets, i.e. fi — 5 x 10 ll) Pa. 
These models of planets are used as starting points to compare 
the influence the ocean depth has on core deformation for dif- 
ferent types of planets. Since we do not try to estimate this de- 
formation for realistic planets, we will not discuss in this section 
the validity of the values we use for the parameters. 
Fig. 5 shows that for a planet with a shallow fluid shell 
(i.e. when R c > a x R p , where a = {0.840,0.915,0.937} for an 
Earth-, Jupiter- and Saturn-like planet respectively), F is less 
than unity, which means that the ocean exerts a loading effect 
on the solid core which is stronger than the gravitational forces 
and opposed to it. That is the case of the Earth where the depth 
of the oceanic envelope does not exceeds 1 % of the size of the 
planet, but giant planets are supposed to have a solid core no 
bigger than the third of the planet size. According to Fig. 5, F 
may reach values up to 2.3 for this kind of planets. That means 
that for a planet with a deep fluid envelope, the ocean tide has 
no loading effect on the core but exerts a gravitational force that 
amplifies the tidal deformation. We refer the reader to Dermott 
(1979) for a complete discussion. 

3.5. Modified Love numbers 

From Eq. (56) we deduce the Love number (cf. Eq. 28) which 
measures the surface deformation: 



(60) 



Quantity 


Value 


Mf (kg) 


5.9736 x 10 24 


KJ (m) 


6371 x 10 3 


Mf (kg) 


Mf - 1.4 x 10 2 


Rf (m) 


Rf, - 3682 x 10 



Reference 



http://nssdc.gsfc.nasa.gov/ 
planetary/factsheet/earthfact.html 



Charette& Smith (2010) 



Table 4. Earth parameters. 



Quantity 


Jupiter 


Saturn 


Reference 


R„ (m) 


10.973 x Rf 


9.140 xi?® 


http://nssdc.gsfc.nasa. 


M„ (kg) 


317.830 xMf 


95.159 xM® 


gov/planetary/factsheet/ 



Table 5. Mass and mean radius of Jupiter and Saturn. 





R c (m) 


M c (kg) 


Reference 


Jupiter 


0.126 xR p 


6.41 xA/^ 


http://www.oca.eu/guillot/ 
jupsat.html (Guillot, 1998) 


Saturn 


0.219 xR p 


18.65 x Mf 


Hubbard et al. (2009) 



Table 6. Mass and mean radius of Jupiter's and Saturn's cores. 



Earth-like planet 
Jupiter-like planet 
Saturn-like planet 
oeean-free planet 




1.83 



0.965 

0.4 0.6 

RJRp (at fixed R c ) 

Fig. 5. Factor F, accounting for the overloading ex- 
erted by the tidally deformed oceanic envelope on the 
solid core, in terms of the ocean depth throw the ratio 
R c /Rp for three types of planet. Parameters are given in 
Tables 4-5-6: for Earth-, Jupiter- and Saturn-like plan- 
ets, we assume respectively that: Rp = {1, 10.97, 9.14} 



(in units of RJ), M p 



[1,317.8,95.16} (in units of Mf,), 



R c = {0.99, 0.126, 0.219} xR p , M c = {0.33,6.41, 18.65} (in 
units of Mf ) and /i = 5 x 10 10 (Pa) for all cores. Note that this 
figure is similar to that drawned by Dermott (1979), divergences 
coming from the use of different values of the parameters. 



Let us give here the expression of the second-order Love number 
(20) associated with the solid core. First of all, let us recall its 
value for an ocean-free planet: According to Eq. (29) with r = 
R c , we get 

3 3 1 
k2 = -h 2 = -—^. (61) 

5 2 1 +yU 

In presence of an ocean on top of the solid core, we may also 
introduce the modified Love number 



def <D'(/? f ) V'(R C ) - U(R C ) V'(R C ) 



U(R C ) 



U(R C ) 



U(R C ) 



1 



(62) 



where V'(R C ) and U{R C ) are obtained from Eqs. (31b) and (21) 
respectively, with r = R c . Thus, expressing £ C /R C in function 
of the modified second-order Love number according to 
Eq. (28), and using Eq. (45), we obtain the expression of in 
terms of h F ~ : 



Id - 
K 2 ~ 



l + 2 -- 

5 1 



(63) 



In this section we have studied the impact of the presence of a 
fluid envelope in the determination of the deformation imposed 
on an elastic core under tidal forcing. In the following, we con- 
sider that the solid core also presents viscous properties so that 
its response to the tidal force exerced by the perturber is no more 
immediate, thus inducing dissipation. The next section addresses 
the quantification of this conversion of energy, which drives the 
dynamical evolution of the whole system. 

4. Anelastic tidal dissipation: analytical results 

Assuming that the anelasticity is linear, the correspondence prin- 
ciple established by Biot (1954) allows us to extend the formu- 
lation of the adiabatic elastic problem to the resolution of the 
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4.2. Case of a two-layer planet 
Let us introduce the quantities 



equivalent dissipative anelastic problem. For initial conditions 
taken as zero and similar geometries, the Laplace and Fourier 
transforms of the anelastic equations and boundary conditions 
are identical to the elastic equations, if the rheological param- 
eters and radial functions are defined as complex numbers. We 
will then denote 

cr = &\ +i&2 (64) 
the complex stress tensor, and 

e = i\ + ie 2 



the complex strain tensor. 

The perturbative strain is cyclic, with tidal pulsations cr 2 , m ,p,q- 
For sake of clarity, we will use from now on the generic notation 
<w = o~2,m,p,q, recalling that there is a large range of tidal frequen- 
cies for each term of the expansion of the tidal potential. The 
stress and strain tensors take the following form: 



o-{oS) = (& x + i& 2 )e iu " , 

e(w) = ioc' W . 
The complex rigidity 

p(cS) = fn(aS) + i p 2 (cS) ■ 



(66a) 
(66b) 

(67) 



where p\ represents the energy storage and p 2 the energy loss of 
the system, is defined by 



p(u) = 



o-(^) 

i(w) 



We may also define the complex effective rigidity 
p.(u>) = p.\(oS) + ip 2 ((jS) 

by: 



where (see Eq. 27) 



p{aj) = yfi((A)) , 



19 



= ^ = ^ = 

yU jU 2p c g c R c 



(68) 

(69) 
(70) 

(71) 



4.1. Case of a fully-solid planet 

The complex Love number k 2 may be expressed in terms of the 
complex effective rigidity p., by: 



kiipS) = - 



1 



2 1 + fi(cj) 

3 1 



(72) 



2 1 +y [//i(w) + ip 2 (to)] 



in the case of a completely solid planet. 

The real part of k 2 characterizes the purely elastic deforma- 
tion, since —Im (k 2 ) gives the phase lag due to the viscosity. 
Therefore, we define the factor of tidal dissipation Q, that repre- 
sents the dissipation rate due to viscous friction, by 



Im £2(01) 

' \h(oj)\ 



Then, from equation (72), we deduce that 
Q(co) = Jl + 



1 



p 2 (CO) H2(U>) 



(73) 



(74) 



1 + 



Po 
Pc 

Po 3 Po 



2a 



and: 



B= + \\- r —\ — 



2p e 



Po 



Pc 



9 (R 



4a \R 



Pc 



(65) in the expression (56) of F: 



A(l+fi) 

F = , 

B + p 

and let us define its complex equivalent 
- A(l+p) 



B 



(75a) 
(75b) 

(76) 

(77) 



The determination of how the presence of an oceanic envelope 
will modify the tidal dissipation consists in the determination of 
k~ , defined as the complex Love number k 2 in presence of the 
fluid envelope. According to the correspondence principle, this 
number is given by the complex Fourier transform of (63), i.e. : 



11 (to) = 

From (63) we get then 
1 



1 + - 



2 p(a>) 



5 1-* 

Pc > 



hUu)- 1 



(78) 



(B + p\) + p\ 



(B+p 1 )\C+—p 1 



3 _ 2 
+ Ya^ 



iA D fi 2 



(79) 



where we made use of the dimensionless quantities a, A and B 
previously defined (see respectively Eqs. 54 and 75), C and D 
given by: 



2 



1-^ 



p c 2a 4a\R„ 



1 + 



1-^ 



2a \R 



(80a) 
(80b) 



Finally, the dissipation factor Q, defined here by 



(81) 



is of the form: 



1 + 



\ 4a 2 A 2 D 2 



1 + 



(B +fii(io)) {^f+p^oS)) 



p 2 (u>) 



(82) 

Thanks to the correspondence principle, one is able to derive this 
general expression of the tidal dissipation valid for any rheology. 
The obtained formulae explicitly reveal its dependence on the 
tidal frequency a> = a- 2 m p q , as shown, for example, by Remus 
et al. (2012), Ogilvie & Lin (2004-2007) for fluid layers. 
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4.3. Implementation of an anelastic model 

Since the factor Q depends on the real and imaginary compo- 
nents jlx and p.2 of the complex effective shear modulus ft, we 
need to define the rheology of the studied body to express it in 
terms of the constitutive parameters of the material. 
The anelasticity of a material is evaluated by a quality factor Q a 
defined by 



(83) 



We may express the solid tidal dissipation, given by Q (Eq. 74) 
for a fully-solid planet and Q (Eq. 82) in the case of a two-layer 
planet, in terms of this factor: 



Q(to = o- 2 , m , p , q ) = \\ 1 



1 



y"2(w) 



Qa(0>) 



(84) 



and 



Q(cj = a-2, m , p , q ) 



II + 



4a 2 A 2 D 2 



1 + 



3yU 2 (w) 



Ga(w) 



(85) 



All previous results are independent of the viscoelastic Theolog- 
ical model. We have now to apply these general expressions to 
a specific rheology which will depend on the physical properties 
of the material. Considering our lack of knowledge on the in- 
ternal structure of giant planets, we will implement the simplest 
model, namely the Maxwell model, which presents the advan- 
tage to involve only two parameters and thus to be easy to use 
(Tobie 2003, and Tobie et al. 2005). A critical overview of the 
four main rheological models has also been done by Henning et 
al. (2009), thus we refer the reader to the three above mentioned 
papers for a detailed comparison. 

4.4. The Maxwell model. 

This model considers a viscoelastic material as a spring-dashpot 
serie. The instantaneous elastic response is characterized by a 
shear modulus G , and the viscous yielding is represented by a 
viscous scalar modulus 77 (see Fig. 6). Notice that the shear mod- 
uli G and fi (introduced in Sect. 2.4) designate the same quan- 
tity. We change here the notation to avoid any confusion with the 
complex shear modulus jl used to study the anelastic tidal dissi- 
pation, and whose real and imaginary parts involve both G and 

n- 




Fig. 6. Representation of the Maxwell model and its correspond- 
ing notations. 

The constitutive equation is given by Henning et al. (2009): 

G (t(cj) + 77 ct(cj) -Grj e(cJ) , (86) 



where the time derivative of a given quantity is denoted by a dot. 
Recalling the relation (68) this equation becomes: 



Ge(u)fi(u)) + 77 — 



fi(<jS) e(o>) 



-Grj e(cj) . 



(87) 



Therefore the real part n\ and the imaginary part fi2 of the com- 
plex shear modulus fi have the following expressions: 



yUi(w) 



yU 2 (w) 



l 2 Go? 



G 2 



T] 2 OJ 2 



77 G oj 



G 2 



rj 2 to 2 



The anelastic quality factor Q a is then given by 

Qa(oJ) = ——; = -pr = ojt m , 
^2(io) G 



(88a) 
(88b) 

(89) 



where tm = rj/G is the characteristic time of relaxation of the 
Maxwell model. As confirmed by Fig. 7, Eq. (89) shows that Q a 
increases linearly with the frequency of the cyclic tidal strain: 
the shorter the oscillation period, the lower the dissipation due 
to intrinsic viscoelastic properties of the material. Moreover, the 
anelastic quality factor is proportional to tm = r//G, so that it 
dissipates more if it is more rigid and less viscous. 



1.X10 3 




4.xl(T 4 6.x 10-" 
oj (rad-s~') 



— 1; = to 3 Pa-s 
-11= 10" Pa-s 
-- i] = 10" Pa-s 
8.x lO" 4 1.x lO" 3 



Fig. 7. Anelastic quality factor Q a of the Maxwell model in 
function of the tidal pulsation a> for different values of the viscos- 
ity 77. G is taken equal to 5 x 10 Pa (see Henning et. al 2009). 
Q a is represented on a logarithmic scale. 

Thus, we may express 112 (Eq. 88b) in terms of the anelastic qual- 
ity factor Q a (Eq. 89): 



J"2(w) = 



(CJT M ) +L0T M 



(90) 



In the case of a fully-solid body, we get, from Eqs. (83) and (90), 
the imaginary part of the complex Love number k.2 (72): 



t \1 i ^ 3r Gojtm 

Im \k 2 (oj)\ = -— -J 

J 2 I + (■ ■ 



(91) 



(ojt m ) 2 (1+yG) 2 
Therefore the dissipation factor Q defined by (74) is of the form: 

Q(a> = 0-2,m,p,q) = 

-^1 + |^ [(uT M y l + wt m ] + wt m | . (92) 
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In the more general case of a two-layer body, the imaginary part 
of the complex Love number given by Eq. (79), takes a dif- 
ferent form than Eq. (91) because of the presence of the fluid 
envelope, so does the two-layer dissipation factor Q (82) with 
respect to its oceanless form (92). To obtain them, one needs to 
replace the shear modulus jj and the anelastic quality factor Q a 
by their expression in the case of the Maxwell model (Eqs 88 
and 89 respectively). 

Fig. 8 compares the dissipation of the solid core with and without 
the presence of a fluid envelope of variable depth for a Saturn- 
like planet, using the parameters given by Tables 5-6-7: as ex- 
pected, the difference between the two decreases with the size of 
the fluid envelope down to about 0.34 x R p ; but for a thiner fluid 
shell, the dissipation get lower than it would be without it. 



0.6 
0.5 
0.4 
$^ 0.3 
0.2 
0.1 



0.0 



X 1 






MM 
















CM 










r t ' 1 



0.2 



0.4 0.6 

RJR P (at fixed R c 



Fig. 8. Relative difference between an d |^2|/2- We 

may distinguish two regimes: for < R c < 0.661/?;, 
the presence of the fluid envelope increases the tidal dis- 
sipation; over this value, tidal dissipation is lower than 
it would have been without fluid shell. Parameters are 
given in Tables 5-6-7 for a Saturn-like planet perturbed at 
the tidal frequency of Enceladus (to — 2.25 x 10~ 4 rad ■ s _1 ): 



R P =9.UF%, M p 



95.16 M® R r 



G = 5 x 10 10 (Pa) and 77 = 10 15 (Pa ■ s). 



0.219 R®, M c 

p 



18.65 M®, 



Planet 


Jupiter 


Saturn 


Satellite 


Io 


Enceladus 


to (rad • s" 1 ) 


2.79 x irr 4 


2.25 x irr 4 


Reference 


Ioannou et Lindzen (1993) 


Lainey et al. (2012) 



Table 7. Tidal frequencies considered in numerical applications. 



present knowledge of extrasolar giant planets, and also planets 
of our Solar System like Jupiter or Saturn, suffers some uncer- 
tainties on the values of these parameters, so that the range of 
values taken by core properties of giant planets presents poor 
constraints. Moreover, even if the presence of a core in Jupiter is 
not yet confirmed (see Guillot 1999-2005), but new data coming 
from seismology may provide more constraints on giant plan- 
ets internal structure (see Gaulme et al. 201 1). Nevertheless, we 
will explore the resulting tidal dissipation of such bodies around 
values of the structural and rheological parameters taken as ref- 
erence and corresponding to those of the literature. 

5. 1 . Baseline structural and rheological parameters 

As reference models, we chose Jupiter and Saturn although their 
core parameters are still uncertain. The values of the global sizes 
and masses of these planets are those of Table 5. 



5.1 .1 . Size and mass of the core 

Presently, two main types of models are available for Jupiter's 
interior. The NHKFRB group 1 uses a three-layer model with 
a thin radiative zone, close to previous models by Saumon 
& Guillot (2004), whereas the MHVTB group 2 proposes a 
new type of Jupiter model that possesses only two layers 
(see Militzer et al. 2008). But, as explained in Militzer & 
Hubbard (2009), the crucial difference lies in the treatment 
of the molecular-to-metallic transition in dense fluid hydro- 
gen, that leads to very different conclusions. The first group 
predicts a small core of less than 10 M® (Saumon & Guillot 
2004), while the second obtains a larger core of 14-18 M p 
(Militzer et al. 2008). Among all these models of Jupiter's in- 
terior, we choose as reference the adiabatic model with Plasma 
Phase Transition (PPT) (http://www.oca.eu/guillot/jupsat.html) 3 
of Guillot (1998), which is of the first type. It predicts a core 
of radius R c = 0.126 x R p and mass M c = 6.41 x M®. Only the 
mass of the core of this reference model will be used in what 
follows. The core radius will serve just as a first approximation, 
as a starting point in our study, since we will present our results 
for several core sizes. 

There are also different models of Saturn's inte- 
rior. According to the model of Guillot with PPT 
(http://www.oca.eu/guillot/jupsat.html) 3 , Saturn's core may 
have a mass of M c = 6.55 x Mf } and a size of R c - 0.174 x R p . 
More recently, Hubbard et al. (2009) infer, from Cassini- 
Huygens data, that Saturn has a larger core in the range 
M c = 15-20 x M p and a corresponding radius of more than 20% 
of the planet size. We adopt this latter as reference model of 
Saturn, with M c = 18.65 x M e and R c = 0.219 x R p . 
All these models support core accretion as the standard process 
for the formation of giant planets. The corresponding parameters 
are listed in Table 6. 



5. Anelastic tidal dissipation: role of the structural 
and rheological parameters 

With our choice of the Maxwell model to represent the rheol- 
ogy of the solid parts of the planet, the dissipation quality factor 
Q depends on the tidal frequency u> and on four structural and 
rheological parameters: the relative size of the core (R c /R p ), the 
relative density of the envelope with respect to the core (p„/p c ), 
the shear modulus (G) and the viscosity of the core (rj). Our 



5.1 .2. Rheological parameters of the core 

The main uncertainties concern the viscoelastic properties of 
the core, namely its shear modulus G and its viscosity 77. At 

1 Nettelmann, Hoist, Kietzmann, French, Redmer and Blaschke 
(Nettelmann et al. 2008). 

2 Militzer, Hubbard, Vorberger, Tamblyn, and Bonev (Militzer et al. 
2008). 

3 It is constructed with CEPAM, Code d'Evolution Planetaire 
Adaptatif et Modulaire (Guillot & Morel, 1995). 
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TIDAL EFFECTIVE DISSIPATION FACTOR Q D£i OF A JUPITER-LIKE PLANET (Mj, = 6.41 X M* ) AT THE TIDAL FREQUENCY OF 10 - 2,79x10 4 rad ■ B 1 ) 
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Fig. 9. Dissipation quality factor Qeff of the Maxwell model in function of the viscoelastic parameters G and 77. Top: for a Jupiter- 
like two-layer planet tidally perturbed at the Io's frequency a> - 2.79 x 10~ 4 rad ■ s . Bottom: for a Saturn-like two-layer planet 
tidally perturbed at the Enceladus' frequency to 2.25 x 10 4 rad ■ s . The red dashed line corresponds to the value of Q e ff = 
{(3.56 ± 0.56) x 10 4 , (1.682 + 0.540) x 10 3 } (for Jupiter and Saturn respectively) determined by Lainey et al. (2009-2012). The 
blue lines corresponds to the lower and upper limits of the reference values taken by the viscoelastic parameters G and rj for an 
unknown mix of ice and silicates. We assume the values of R p = {10.97, 9.14} (in units of i?®), M p = {317.8, 95.16} (in units of M p ), 
M c = {6.41, 18.65} x given in Tables 5-6. 



high pressure and temperature, theoretical models and experi- 
ments show that G and 77 values depend on temperature and pres- 
sure. But no experiments are available at the very-high pressures 
and temperatures we may expect in Jupiter's and Saturn's cores 
(Guillot 2005). Nevertheless, geophysical and experimental data 
have allowed to constrain the rheology of the icy satellites of 
Jupiter, since their ranges of pressure and temperature are simi- 
lar to those of the outer mantle of the Earth (Tobie 2003). Then, 
keeping in mind that these values may differ by several orders of 
magnitude in our case, we will adopt reference values based on 
these data, assuming that Jupiter's and Saturn's cores are made 
of ice and rock. We will then explore, in the following Sect. 5.2, 
the variation of the tidal dissipation for a large range of values 
of the rheological parameters around those taken as reference. 
We thus assume that the shear modulus G is in the range 
[G ice = 4 x 10 9 (Pa) , Gsuicau = 10 11 (Pa)] (Henning et al. 2009). 
Concerning the viscosity r\, it takes its values in the range 
T]i ce = 10 14 Pa ■ s,Ti s mcate = 10 21 Pa ■ sj for the icy satellites of 
Jupiter at high pressure (Tobie 2003). We expand this range, re- 
ducing its lower boundary by two orders of magnitude, following 
the discussion of Karato (201 1) which seems to indicate that, at 
the very high pressures, viscosity in deep interior of super-Earths 
may decrease by two or three orders of magnitude. We refer 



the reader to the Karato's paper for an overview of all plausi- 
ble mechanisms that may cause a change in the viscous-pressure 
relationship at very-high pressures. 



5.2. Dependence of tidal dissipation on rheology 

Since tidal dissipation causes exchange of angular momentum 
in the system, it may be quantified by monitoring carefuly the 
orbital motion of the system. Using astrometric data covering 
more than a century, Lainey et al. (2009-2012) succeeded in 
determining from observations the tidal dissipation in Jupiter 
and Saturn: namely, Qjupiter = (3.56 + 0.56) x 10 4 determined by 
Lainey et al. (2009), and g S atum = (1-682 + 0.540) x 10 3 deter- 
mined by Lainey et al. (2012) and requireded by the formation 
scenario of Charnoz et al. (201 1). However, with such a method, 
the different contributions to the global tidal dissipation, com- 
ing from each layer constituting the planet, are lumped together. 
Equations of the dynamical evolution (Eqs. 99 to 103) link the 
observed evolution rates of the rotational and orbital parame- 
ters to the observed tidal dissipation and system characteristics. 
Since all these rates are proportional to R p , where R p is the planet 
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oj (rad.s _1 ) RJRp (at fixed R p ) 



Fig. 10. Dissipation quality factor £) e ff normalized to the size of the planet for Jupiter-like and Saturn-like giant planets. Note that all 
curves are represented with a logarithmic scale. Left: dependence to the perturbative strain pulsation to, with R c = {0.20, 0.34} x R p . 
Right: dependence to the size of the core, with at 2.25 x 10~ 4 rad ■ s _1 (tidal frequency of Enceladus) for the blue curve associated 
to a Saturn-like planet, and a> 2.79 x 10~ 4 rad ■ s _1 (tidal frequency of Io) for the red curve associated to a Jupiter-like planet. The 
green curve corresponds to the prescription of Goodman & Lackner (2009) (see plain text for details). The red and blue dashed 
lines correspond to the value of Qeff = {(3.56 + 0.56) x 10 4 , (1.682 + 0.540) x 10 3 } (for Jupiter and Saturn respectively) determined 
by Lainey et al. (2009-2012). Their zone of uncertainty is also represented in the corresponding color. We assume the values of 
R p = {10.97, 9.14} (in units of R$), M p = {317.8, 95.16} (in units of M®), M c = {6.41, 18.65} x M® given in Tables 5-6. We also take 
for the viscoelastic parameters G = {4.85,4.45} x 10 10 (Pa), and r\ = {1.26, 1.78} x 10 14 (Pa • s _1 ) for Jupiter and Saturn respectively. 



radius, we introduce the associated dissipation factor 4 : 



H(R P 
x — 



m(R e )\ 



(93) 



where \k£(R c )\ can be deduced from Eq. (81), and \k^(R p )\ des- 
ignates the modulus of the second order Love number of the 
planet's surface that is obtained from Eqs. (52) and (21) 



E(R P ) 



V'(R P ) 
U{R P ) 



1 



ft) 



2 Po 



where H acounts for the quantity 



(94) 



Therefore, we can evaluate the real part of the second 
order Love numbers k^iRc) and k^{R p ), accounting re- 
spectively for the deformation of the core's and planet's 
surface, for parameters whose values are compatible with 
the tidal dissipation observations (Lainey et al., 2009- 
2012). For Jupiter and Saturn, in this order, assuming 
that R c = {0.170, 0.260} XR P , G = {4.85,4.45} x 10 10 Pa 
and r] = {1.26, 1.78} x 10 14 Pa • s, we obtain that 
'Re [k%(R c )\ = {3.21,3.31} and 'Re [k^(R p )] = {1.37,0.24}. 
These estimations at the planet's surface can be compared to the 
value of Gavrilov & Zharkov (1977) of ko = 0.379 for Jupiter 
and k2 = 0.341 for Saturn obtained for stratified models. As 
discussed in the aforementioned paper, the differences between 
both evaluations are linked to the degree of stratification: the 
more the planet interior is stratified, the smaller the second order 
Love number (we recall that the second order Love number of a 
homogeneous fluid planet is 3/2). 



Po 
Pc 



3po 
2p e 



Po 
Pc 



3 R 
2\R 



Po 
Pc 



(95) 



Since we have weak constraints on the viscoelastic parameters of 
giant planet cores (Guillot 2005), we thus have to explore a large 
range of values. Fig. 9 shows the tidal dissipation factor Q e ff 
around the reference values presented in Sect. 5.1, expanding 
the range up to about + 2-4 orders of magnitude for G and rj. In 
the middle region (inside the blue rectangle on Fig. 9), where tj 
and G correspond to the reference values, the dissipation factor 
2etf of Saturn (resp. Jupiter) may reach values of the order of 
10 3 (resp. 10 4 ), and in the whole field it varies up to 10 20 . 
From Fig. 9, we deduce that the tidal dissipation of the core may 
reach the values observed for Jupiter (Lainey et al. 2009) and 
Saturn (Lainey et al. 2012) assuming that Jupiter's core (resp. 
Saturn's core) has a radius 34.92% (resp. 18.72%) larger than 
this of Guillot 1998 (resp. Hubbard 2009). 



4 This can also be demonstrated by calculating the perturbation of 
the gravific potential at the surface of the planet, adapting Zahn (1966) 
theory taking into account the boundary conditions at the core surface. 



5.3. Dependence of tidal dissipation on the size of the core 
and the tidal frequency 

For Fig. 10, the values of the viscoelastic parameters G and rj, 
and of the core size R c result from Fig. 9: we chose them so as the 
tidal dissipation factor Q reaches the observed values of Lainey 
et al. (2009-2012) on the condition that the rheological parame- 
ters stand in the more realistic domain defined by the lowest and 
highest value of the ice and rock viscoelasticities taken as refer- 
ence. Taking into account the global dissipation values obtained 
by Lainey et al. 2009 & 2012 for Jupiter and Saturn, we may 
deduce some constraints on the viscoelastic parameters and also 
the size of the core (looking at the red dashed line in Fig. 9). We 
thus will assume that they take values allowing such a dissipa- 
tion: 

- we first chose a core slightly larger than that assumed until 
now: R c = 0.170 x R p for Jupiter, and R c = 0.260 x R p for 
Saturn, 

- we then fixed the value of the shear modulus G to the low- 
est value needed to reach the observed tidal dissipation of 
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Lainey etal. (2009-2012), i.e. G = 4.85 x 10 10 Pa for Jupiter, 
and G = 4.45 x 10 10 Pa for Saturn, 
- e finally searched the more realistic value of the viscos- 
ity which corresponds to the observed tidal dissipations of 
Jupiter and Saturn: tj — 1.26 x 10 14 Pa ■ s for Jupiter, and 
j] = 1.78 x 10 14 Pa • s for Saturn. 

Fig. 10 explores, for the present model, the dependence of Q e ff 
on the pulsation u> and the size of the core R c normalized by the 
size of the planet R p . With these parameters, the figure indicates 
that Saturn dissipates slightly more than ten times more than 
Jupiter, since (/? c ) Satum > (fl c )j up iter and (p c ) Satum < (p c ) ju pi, el - In 
the range of tidal frequencies of Jupiter's and Saturn's satellites 
(2.25 x 10~ 4 rad ■ s" 1 < u < 2.95 x 10~ 4 rad ■ s"\ Lainey et 
al. 2009-2012), the effective dissipation factor Q^s keeps the 
same order of magnitude. However, it strongly depends on the 
size of the core, as it may loose up to 6 orders of magnitude 
between a coreless planet and a fully-solid one. Note that for 
a given core (so that M c , R c and then p c are fixed) and a given 
mass of the planet M p , the density of the fluid envelope p a , which 
varies with its height R p - R c , can not exceed p c . Since: 



Po(Rc) = 



Mp - M e 
4/3 7r(R 3 p -Rl) ' 



(96) 



this condition gives a limit for the core size: 



P ' sup 



Ma 

M„ 



1/3 



(97) 



In 2004, Ogilvie & Lin also studied tidal dissipation in giant 
planets, and particulary the tidal dissipation resulting from the 
excitation of inertial waves in the convective region by the tidal 
potential for rotating giant planets with an elastic solid core. 
They obtained a decrease of the quality factor Q of one or- 
der of magnitude considering the dynamical tide (due to iner- 
tial modes) compared to the equilibrium one: from Q — 10 6 to 
Q = 10 5 , but it is not efficient enough to explain the observed 
tidal dissipations in Jupiter or Saturn which are of 1-2 orders 
of magnitude higher (Lainey et al. 2009-2012). Moreover, they 
showed that the dissipation resulting from the resonance be- 
tween fluid tide and inertial modes highly depends on the tidal 
frequency in the range of inertial waves, as do also the core- 
less models (Wu 2005). This disagrees with the weak frequency- 
dependence obtained from astro me try (Lainey et al. 2012). 
By discussing the size of the core, Goodman & Lackner (2009) 
got a higher value of the quality factor Q, in the range 10 7 - 

10 8 x (0.2 R p /R c ) which disagrees with the observed value of 
the tidal dissipation of Saturn (see Fig. 9). 
The present two-layer model proposes an alternative process 
to reach such a high dissipation with a smooth frequency- 
dependence of Q. But the uncertainties on the structural and The- 
ological parameters does not allow us to firmly conclude that the 
tidal dissipation of the core can explain by its own the tidal dis- 
sipation observed in giant planets of our Solar System by Lainey 
et al. (2009-2012). 

Through these expressions of the tidal dissipation closely linked 
to the internal structure of the planet and its rheological proper- 
ties, we are now able to derive the equations of the dynamical 
evolution of the system with explicit dependence on the tidal 
frequency. 



5.4. Comparison with previous work of Dermott 

The difference beetwen our study and the work of Dermott 
(1979) lies in the treatment of the tidal dissipation. 
Dermott draws his expression from an evolution scenario of 
Satumian and Jovian systems. He assumes that the satellites of 
Jupiter and Saturn were formed 4.5 x 10 9 ago, and that their 
semi-major axis have changed by 10 % since their formation, 
with a stable resonance of the main satellites of Jupiter (Io, 
Europa and Ganymede) since their formation but a young res- 
onance of the satellites Mimas and Thetys of Saturn. He also 
assumes an average value of the tidal dissipation, independent 
of amplitude, frequency and time. All these assumptions lead 
Dermott (1979) to a tidal dissipation factor Q that only depends 
on the mass M c , the size R c , the elasticity p. of the core and a 
dimensionless coefficient K characteristic of the evolution sce- 
nario of the planet. In particular, his tidal factor Q is directly 
proportional to R s c (Eq. 27 of Dermott 1979), so that the tidal 
dissipation gets lower as the core size increases (see Fig. 4 of 
Dermott 1979), while one should expect an opposite behaviour. 
Instead, our model is constructed on physical considerations of 
the internal structure and properties of the core. In particular, we 
derived our tidal dissipation factor Q with no assumption on the 
evolution of the Jovian and Satumian systems. To do so, we used 
the correspondence principle of Biot (1954). This allowed us to 
obtain an expression of the tidal dissipation factor valid for any 
rheological model of planets' cores. Moreover, our expression 
(82) of Q not only depends on the mass M c , the size R c and the 
elasticity p of the core, but also on the tidal frequency u> and 
the viscosity r\. We notice, in particular, that tidal dissipation in- 
creases when the size of the core increases, contrary to Dermott's 
strange result. 

6. Equations of the dynamical evolution 

Mass redistribution due to the tide generates a tidal torque of 
non-zero average which induces an exchange of angular momen- 
tum between the orbital motion and the rotation of each com- 
ponent. As shown in MLP09 & Remus et al. (2012), this tidal 
torque is proportional to the tidal dissipation ratio )M \/Q (see 
also Correia & Laskar 2003, Correia et al. 2003, and Murray 
& Dermott 2000). Note that for a perfectly elastic material, the 
core will be elongated in the direction of the line of centers, in- 
ducing a torque, I r A (p f VL0 d'V, with periodic variations of 

zero average, so that no secular exchanges of angular momentum 
will be possible (see Zahn 1966a and Remus et al. 2012). On the 
other hand, if the core is anelastic, the deformation of the core 
resulting from the equilibrium adjustment presents a time delay 
At with respect to the tidal forcing, which may be measured also 
by the tidal lag angle 26/ or equivalently by the quality factor Q 
(see Ferraz-Mello et al. 2008 or Efroimsky & Williams 2009): 



tan [A? x cr 2 , m , M ] = tan \26 (cr 2 , miM )] 



1 



Q (f 2,m,p,q) 



(98) 



Thus, the tidal bulge is no more aligned with the line of centers, 
as shown in Fig. 11. 

The resulting tidal angle is at the origin of a torque of non-zero 
average which causes exchange of spin and orbital angular mo- 
mentum between the components of the system. 
The evolution of the semi-major axis a, of the eccentricity e, of 
the inclination /, of the obliquity s and of the angular veloc- 
ity CI (Ia denotes the moment of inertia of A), is governed by 
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where the functions 

d 2 



m ,j lPl q(e,I,e) are expressed in terms of 
;m ( e X Fi,j,p(J), and G2, P ,q(e), which are defined in Sect. 2.2 



< Hxm,j, P ,qi.e, I, e) = 



5 (2-1/1)1 
Un(2 + \j\)\ 



xd 2 jm (e)F 2J , p (I)G 2 , M (e), (104) 



Fig. 11. Tidal interaction involving a solid body. Body B exerts 
a tidal force on body A, which adjusts itself with a phase lag 
2 6, because of internal friction in the anelastic core. This ad- 
justment may be split in an adiabatic component, corresponding 
to the elastic deformation, in phase with the tide, and a dissipa- 
tive one, resulting from the viscous internal frictions, which is in 
quadrature. 



the following equations, established in MLP09 and Remus et al. 
(2012): 



and /? eq designates the equatorial radius of body A. 

From these equations one may derive the characteristic times of 

synchronization, circularization and spin alignment: 
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7. Conclusion 

To conclude, the purpose of this work was to study the tidal dis- 
sipation in a two-layer planet consisting in a rocky/icy core and 
a fluid envelope, as one expects to be the case in Jupiter, Saturn, 
and many extrasolar planets. We considered the most general 
configuration where the perturber (star or satellite) is moving on 
an elliptical and inclined orbit around the planet which rotates on 
an inclined axis. We expanded the tidal displacement in Fourier 
series and spherical harmonics, each term of the expansion hav- 
ing a radial part proportional to that of the corresponding term 
of the tidal potential, which depends on the eccentricity, the in- 
clination and the obliquity. We followed the method by Dermott 
(1979) to derive the modified Love numbers /jf and account- 
ing for the tidal deformation at the boundary of the solid core. 
Like Dermott, we made the simplifying assumption that core and 
envelope have a constant density. Then, generalizing the results 
of his work invoking the correspondence principle, we obtained 
the tidal dissipation rate of the core expressed by \kS \/Q, Q be- 
ing the quality factor. That ratio depends on the tidal frequency 
and on the rheological properties of the core; unlike Dermott we 
made no assumption on the formation history of the system. As 
mentioned in Sect. 5.1 the rheological properties of planetary 
cores are still quite uncertain. However, taking plausible values 
for the viscoelastic parameters G and 77, we obtain a tidal dissi- 
pation which may be much higher than for a fully fluid planet 
and weakly frequency-dependent. Under these assumptions, we 
find that the low value of Q = (1.682 + 0.540) x 10 3 , determined 
by Lainey et al. (2012) and needed by Charnoz et al. (2011) to 
explain the formation of all mid-size moons of Saturn from the 
rings, can be reached taking into account the tidal dissipation of 
Saturn's core. In the same way, the dissipation in Jupiter's core 
may explain the value of the Q-factor determined by Lainey et 
al. (2009), i.e. Q = (3.56 + 0.56) x 10 4 . But to do so, we need to 
assume a core in Jupiter and Saturn slightly larger than those of 
Guillot (1998) and Hubbard et al. (2009). But we recall that in 
our model the density was assumed to be piecewise constant. In 
the future, we shall consider a non constant density profile, to 
2 ) evaluate the impact on our results of a realistic density stratifi- 

[Hz, mdp ,q{e,I,s)\ , (103) cation. Moreover, there are much uncertainties on the determi- 
J nation of the core size in giant planets as Jupiter and Saturn, so 



+ q 
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that we need more constraints on the system formation by core 
accretion (Pollack et al. 1996) and differenciation resulting from 
the internal structure evolution (Nettelmann 201 1). Furthermore, 
seismology seems to offer an interesting way to improve our 
knowledge of giant planets interiors (Gaulme et al. 201 1). 
To conclude, the purpose of this paper was to study tidal dissi- 
pation in the solid parts of a simple model of two-layer planet; 
we show how this mechanism may be powerful. The results de- 
rived here are general in the sense that no specific rheological 
model has been assumed. However, considering the lack of con- 
straints on the rheology of giant planets cores, we have chosen 
the simplest Maxwell model to illustrate the tidal dissipation. 
This work represents thus a first step for further numerical inves- 
tigations in more realistic cases. 
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In Sect. 2.3 we gave the system of equations (1 1), with its bound- 
ary conditions (12), governing an elastic planet under tidal per- 
turbation. Considering that such perturbations are of small order 
of magnitude compared to the hydrostatic equilibrium, we pro- 
posed in Sect. 2.4 a method to linearize the system (14). Thus, 
assuming the expansion (18) of all quantities in spherical har- 
monics, we obtain the following system of equations governing 
the scalar radial parts of these expansions (Alterman et al. 1959; 
Takeuchi & Saito 1972): 
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Solutions of (109) are given by (19). 
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